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ABSTRACT 



^ . We discuss 7-ray emissions from the outer gap accelerators of middle-aged 

pulsars for part of the series of our studies. A two-dimensional electrody- 
namic model is used to solve the distribution of accelerating electric fields with 
, electron and positron pair creation and radiation processes in the magnetic 

. meridional plane. We compute the curvature radiation and the synchrotron 

■ 

radiation by solving the evolution of the Lorentz factor and the pitch angle. 
The calculated spectra are compared with observed phase-averaged spectra. 

00 : 

, We also use a three-dimensional geometrical model to discuss the pulse pro- 

files. We argue that the outer gap of middle-aged pulsars occupies the whole 
• 1—1 . 

■ region between the last-open field lines and the critical magnetic field lines, 

H ' 

' which are perpendicular to the rotational axis at the light cylinder. We as- 

sume that there is no outer gap accelerator inside the light cylinder between 
the rotational axis and the critical magnetic field lines. For the Geminga pul- 
sar, we demonstrate that the outward curvature radiation dominates in the 
spectrum above 10 MeV, while the inward synchrotron radiation dominates 
below 10 MeV. We find that the computed spectrum is consistent with the 
observations in X-ray through 7-ray bands. With the pulse morphology of the 
7-ray emissions, we argue that the inclination angle and the viewing angle for 
the Geminga pulsar are a ^ 50° and ^ ~ 90°, respectively. We also apply our 
method to another four middle-aged radio pulsars, whose spin-down power 
and distance from the Earth expect the possibility of detection of 7-ray emis- 
sions from those pulsars by Fermi telescope. Applying the inclination angle 
and the viewing angle inferred from radio polarization characteristics, the vis- 
ibility of the 7-ray emissions from the outer gap is discussed. We show that 
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7-ray emissions from PSRs B0355+54, B1449-64 and B1929+10 will probably 
be detectable by Fermi telescope. For PSR B0740-28, on the other hand, 7-ray 
emissions will not be detected because the 7-ray beam from the outer gap will 
be oriented in a different direction from the viewing angle. 

Key words: pulsars: middle-aged pulsars- radiation mechanisms: non-thermal 

1 INTRODUCTION 

The Energetic Gamma-ray Experiment Telescope (EGRET) aboard the Compton Gamma- 
ray observatory detected the gamma-ray emissions from six pulsars (Thompson 2004), five 
of which show pulsed radio emissions. The sixth, Geminga pulsar, has known as a radio quiet 
7-ray pulsar or with very weak radio emissions (Vats et al. 1999). Because of the observed 
high brightness temperature of radio emissions from pulsars, the radio radiation is a coherent 
mechanism. On the other hand, the emission mechanisms in optical through 7-ray bands 
are non-coherent and non-thermal processes, which have been discussed with the polar cap 
model, the slot gap model and the outer gap model in the literature. 

The polar cap accelerator has been considered as the origin of radio emission from pulsars 
(Ruderman & Sutherland 1975). The configuration of the dipole magnetic field in inner 
magnetosphere region provides a natural explanation of position-angle curve of polarization 
of the radio pulse. For 7-ray emissions, however, the polar cap model has to invoke both 
a small inclination angle and a small viewing angle measured from the rotational axis to 
reproduce the observed wide separation of the two peaks appeared in the pulse profiles 
(Daugherty & Harding 1996). 

The slot gap model (Arons 1983) is an extension of the polar cap model. The possibility 
of a high-altitude extension of the accelerator near the last-open field lines was proposed. 
The recent slot gap model explains the spectrum in optical to 7-ray bands for the Crab 
pulsar (Muslimov & Harding 2003, 2004; Harding et al. 2008) and the pulse profiles in 7-ray 
bands for the Vela pulsar (Dyks & Rudak 2003). 

The outer gap accelerator model assumes a strong acceleration between the null charge 
surface of the Goldreich- Julian charge density and the light cylinder (Cheng, Ho & Ruderman 
1986a, b). The Goldreich- Julian charge density is described by pcj ~ —VLBz/2tic with Vl 
being the rotation frequency, the magnetic field component projected to the rotation 
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axis, and c the speed of light. The recent studies of the outer gap model also explain the 
spectra and pulse profiles measured by EGRET (Hirotani 2007; Tang et al. 2008; Takata et 
al. 2008). Although numerous efforts have been done to discriminate emission models with 
7-ray data, the origin of the non-thermal emission remains inconclusive. 

In addition to the 7-ray data of EGRET, emission properties in lower energy bands have 
also been provided. For example, Harding et al. (2002) found that the pulse profile of the 
Vela pulsar in RXTE bands shows four or five peaks in a single period instead of the double 
peak structure measured by EGRET. The detail polarization features of optical pulsed 
emissions from the Crab pulsar have been measured by Kanbach et al (2005). INTEGRAL 
also measured the polarization of the 7-ray emissions from the Crab pulsar for the first 
time (Dean et al. 2008). In the future, furthermore, the polarization of X-ray and soft 7-ray 
emissions from the Crab pulsar will probably be able to be measured by ongoing projects such 
as PoGo (Kataoka et al. 2005) and the Nuclear Compton Telescope (Chang et al. 2007). 
Both the Astro -rivelatore Gamma a Immagini LEggero (AGILE) and recently lunched 
Fermi Gamma-Ray Space Telescope (Fermi, hereafter) will measure more detail features 
of the phase-resolved spectra and the pulse profiles in 1 MeV through 100 GeV bands. In 
particular, Fermi telescope is expected to measure pulsed 7-ray emissions from both many 
radio pulsars and from unidentified EGRET sources. A comprehensive theoretical study 
with observations in multi- wavelength bands is desired and will be important to discriminate 
emission models. 

Takata et al. (2007) and Takata & Chang (2007) discussed non-thermal emissions from 
the young pulsars, PSRs B0531+21 (Crab) and B0540-69, with the outer gap model. In 
Takata et al. (2007), they demonstrated that the outer gap model can explain observed 
polarization properties of the optical pulsed emissions from the Crab pulsar. Takata & Chang 
(2007) extended their calculation of the polarization to higher energy band and predicted 
the results of future observations. Also, they produced observed phase-resolved spectra in 
optical through 7-ray bands for the Crab pulsar. They predicted that 7-ray flux on the Earth 
from PSR B0540-69 is about lO^^^erg/cm^s, which may be a target for Fermi observation. 

Takata et al. (2008) discussed the emission process in optical through 7-ray bands for the 
Vela-like young pulsars. Solving the electrodynamics of the outer gap, they first pointed out 
that inward emission by ingoing particles is an important component below 10 MeV bands, 
while only outward emission contributes to the spectrum above 10 MeV as the traditional 
outer gap model has assumed. They produced the observed multi-peak structure of the pulse 
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profile in X-ray bands of tlie Vela pulsar with the inward and the outward emissions. They 
predicted that the pulse profile at 1-10 MeV has four peaks in a single period, which can be 
checked in future observations. 

We had studied the non-thermal emissions from the Crab-like and the Vela-like pulsars. In 
this paper, therefore, we discuss the non-thermal emissions from the outer gap of the middle- 
aged pulsars for part of the series of our studies. Using a two-dimensional electrodynamic 
model, we will solve the accelerating electric field with electron and positron pair creation 
and radiation processes (section [3]). We will also discuss the X-ray and 7-ray pulse profiles 
using a three-dimensional geometrical model (section Hj). We apply the outer gap model 
to the Geminga pulsar (section l5.ll) . The calculated spectrum will be compared with the 
observed phase-averaged spectrum in X-ray through 7-ray bands, (section [5. 1.2p . We discuss 
the inclination angle, which is angle between the rotational axis and the magnetic axis, and 
the viewing angle with the observed properties of the 7-ray pulse profile for the Geminga 
pulsar, (section l5.1.3p . Also by assuming that the radio beam from the Geminga pulsar is 
oriented in a different direction from the line of sight, we will further constrain the magnetic 
inclination angle (section l5.2p . 

We will also apply the model to another middle-aged radio pulsars (section [5l3|) . whose 
timing properties are similar to the Geminga pulsar. We will show candidates for detections 
by Fermi telescope. With its high sensitivity, Fermi telescope will also be able to identify 
7-ray quiet radio pulsars, although their 7-ray emissions should be detectable in the sense 
of magnitude of the spin-down luminosity and its distance. With the present outer gap 
study, we will argue that PSR B0740-28 is a candidate of the 7-ray quiet pulsars due to 
misalignment of the directions of the 7-ray beam and the line of sight. 



2 MAGNETOSPHERES OF MIDDLE-AGED PULSARS 

Figure [1] shows a schematic picture of the magnetosphere model for middle-aged pulsars. In 
general, we expect that the outer gap can arise between the last-open field lines and the 
magnetic field lines (so called critical field lines, hereafter) on which the Goldreich- Julian 
charge density vanishes at the light cylinder. Between the rotational axis and the critical 
field lines, the null charge surface is located beyond the light cylinder so that the Goldreich- 
Julian charge density does not change its sign all the way from the stellar surface to the 
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Figure 1. Schematic figure of the magnotosphere of middle-aged pulsars. The outer gap (filled region) occupies region between 
the last-open field line and the critical field line The radio wave originates from the inner region of the magnetosphere. The 
optical to 7-ray photons are from the outer gap region. The critical field line (thick solid line) is defined by the magnetic field 
line on which the Goldroich-Julian density vanishes at the light cylinder. 



light cylinder. Therefore, we expect that there is no outer gap in the region between the 
rotational axis and the critical field lines. 

For a static dipole field, the polar angle of the critical field line on the stellar surface in 
the magnetic meridional plane is 
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where a is the magnetic inclination angle, i?* is the stellar radius, Ric = c/Q is the light 
radius, and ^„ is the angle of the null charge surface of the Goldreich- Julian density; 
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We define the trans-field thickness of the maximally extending gap in the magnetic merid- 
ional plane with SO^ax = ^pc — ^c, where 9pc is the polar cap angle. For the static dipole field, 
the polar cap angle is expressed by 



6pc,o = a sin ^ 



sm[(iic 



a] 



J^lc 



(3) 



where 



9^ 



Ic 



tan 



-3 - VgTstM? 



a 



4 tan a J 

For the middle-aged pulsars, we can expect that the thickness of the outer gap becomes 
to be comparable with the maximum thickness SO^ax as follows. The potential drop on the 
stellar surface between the critical magnetic field line and the polar cap rim is 
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Va ~ 1.3 X 10l=^P-3/2pl/2^^ 




) Volt, 



(4) 



which becomes Va ~ 1.7 x 10^^ Volt for the Geminga pulsar. 

We can show that almost whole potential drop Va is required to explain the GeV emissions 
from the Geminga pulsar. For the Geminga pulsar, a spectral cut off around 2 GeV was 
measured (Fierro et al. 1998). We estimate whole potential drop in the gap, Vgap, from the 
cut off energy. We assume that 7-ray photons are emitted via the curvature radiation, and 
that the particles that emit the cut off photons are accelerated by the whole potential drop 
along the magnetic field lines in the gap. The typical strength of the accelerating electric 
field is estimated from E\\ ~ Vgap/Ric- We assume the force balance between the acceleration 
and the radiation back reaction force to calculate a typical Lorentz factor of the accelerated 
particles; F ~ {3RlE^^/2eY^^ ~ {3RlVgap/2eRicY^'^ , where Rc is the curvature radius of the 
magnetic field line. Using the typical energy of the curvature photons E^. ~ ShcT^ /AnRc, the 
estimated potential drop in the outer gap is 



where we used Rc = Ric- We find that the estimated magnitude of the potential drop is 
comparable with the available potential drop defined by equation of (jl]) for the Geminga 
pulsar. The observed 7-ray emissions from the Geminga pulsar imply that most of the region 
between the last-open field lines and the critical field lines will be occupied by the outer gap 
for the middle-aged pulsars. 

For younger pulsars, the outer gap accelerator occupies about 10 % of the whole region 
between the last-open field lines and the critical field lines. For example, the available po- 
tential drop defined by equation of (j4]) for the Vela pulsar becomes Va ~ 2.5 x 10^^ Volt. 
The observed cut off energy around 3 GeV indicates the potential drop in the gap is about 
Vgap ~ 10^^ Volt. Therefore, only about 10 % of the available potential drop is required to 
explain the 7-ray emissions from the Vela pulsar. 

From the dynamic point of view, we can also argue that the outer gap of the middle- 
aged pulsars occupies most of the region between the last-open field line and the critical 
field line. The 7-rays emitted by the accelerated particles may convert into electron and 
positron pairs by the pair creation process with a soft X-ray field from the stellar surface. If 
typical energy of the curvature photons have an energy such that the pair creation condition 
E^Ex ^ {rrieC^Y is satisfied, the outer gap can not develop further due to screening effects 
of the pairs. In such a case, the position of the upper boundary of the gap is determined 
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by the pair creation condition. For the Geminga pulsar, however, the pair creation in the 
outer magnetosphere is not efficient. The temperature of the stellar surface is about Tg ~ 
0.5 MK (Kargaltsev et al. 2005), so the threshold energy of the 7-rays for the pair creation 
is ~ {mf.c^)'^ / kTg ~ 6 GeV, where k is the Boltzmann constant. From equation of (Q, 
the potential drop of Vgap ~ 10^^ Volt in the gap is required to accelerate the particles to 
emit 6 GeV photons. However it will be impossible to accelerate the particles because the 
available potential drop in the outer gap ~ 2 x 10^^ Volt estimated from equation of (jl]) 
is smaller than the required potential drop for the pair creation. The particles can not be 
accelerated in the outer gap to emit 6 GeV photons with the curvature radiation process. 
Therefore, we expect that the outer gap of the middle-aged pulsars can develop to occupy 
the region between the last-open field lines and the critical field lines. As we said in the ffist 
paragraph in this section, we do not expect that the outer gap accelerator arises beyond the 
critical field lines, because the Goldreich- Julian charge density does not change its sign all 
the way from the stellar surface to the light cylinder 

We assume that the radio emission originates from the polar cap accelerator, while the 
emissions in optical through 7-ray bands originate from the outer gap accelerator. Electro- 
dynamically speaking, because of the different flow directions of the electric currents in the 
polar cap and in the outer gap, the current conservation forces that both acceleration re- 
gions do not exit on the same magnetic field lines. In this paper, we anticipate that the outer 
gap accelerator is extending from the last-open field line to the critical field lines, and the 
polar cap accelerator extends in a region between the critical field lines and the rotational 
axis. This assumption will be rejected if the radio emissions take place on the last-open field 
line. In fact, there are no observational evidences. For example, the phenomenological study 
with the polarization characteristics for the radio emissions (Dyks et al. 2004) indicates that 
the outer rim of the radio beam is emitted on the magnetic fields located well above the 
last-open field lines. 

3 TWO-DIMENSIONAL ELECTRODYNAMIC OUTER GAP MODEL 

In this section, we briefly summarize our method of the electrodynamic study of the outer 
gap model. A more detail description of our model can be found in the previous paper 
(Takata et al. 2004). We solve the accelerating electric field with the radiation and the pair 
creation processes in the magnetic meridional plane, which includes the rotation axis and 
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the magnetic axis. In the electrodynamic study, we use a static dipole field to solve the 
structure of the outer gap. We consider a magnetized rotator in which an inclination angle 
a between the rotational axis and magnetic axis is smaller than 90°. 

3.1 Electrodynamics 

The stationary electric potential, $nco; for the accelerating field in the observer frame is 
obtained from (Shibata 1995; Mestel 1999) 



where p(r) is the space charge density, and pcjij) is the Goldreich- Julian charge density, 
for which we use pcj = —QBz/27rc. 

For the steady state, the continuity equations of the particles are written as 



where v\\ ~ c is the velocity along the field line, S{t) is the source term due to the pair 
creation process, and and denote the number density of the outgoing (positrons) 
and ingoing (electrons) particles, respectively. In the outer magnetosphere, 7 + X — > e"*" + e~ 
pair creation process contributes to the source term S{r). We simulate the pair creation 
process with a Monte Carlo method, which is explicitly described in Takata et al. (2004, 



To solve the Poisson equation we impose the boundary conditions on the four bound- 
aries, which are called as inner, outer, upper and lower boundaries. The inner and outer 
boundaries are defined by the surfaces on which the accelerating electric field is vanishes, 
that is, E\\ = 0. On the inner, upper and lower boundaries, the accelerating potential is 
equal to zero, that is, $ = 0. On the inner boundary, because both conditions E\\ = and 
$ = are imposed, we can not know the position of the inner boundary in advance. We 
solve the position of the inner boundary with the dynamics. 

The lower boundary is defined by the last-open field line. For the upper boundary, we 
define with a magnetic field line. For the outer boundary, we define it with a curve line, 
which is perpendicular to the magnetic field lines in the magnetic meridian plane. For such 
a case, the upper part of the our gap extends beyond the light cylinder, if we put the outer 
boundary of the lower part of the gap close to (but inside ) the light cylinder. The present 
static Poisson equation of ([6]) will not give the correct solutions for the electric field outside 
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the light cyhnder. As we will see later (Figure [5]), however, most 7-rays are emitted from 
the lower part of the gap. The spectral features will be less affected by the structure of the 
upper part of the gap. We expect that the results do not change much between the present 
approximated and correct treatments. To calculate the spectra, we take into account the 
emissions inside the light cylinder. 

3.2 Particle motion 

To calculate the synchrotron and curvature radiations, we solve an evolution of the particle's 
momentum using the electric field distribution in the outer gap obtained by the method 
described in section [XT! The equation of motion for momenta of the parallel {P\\/meC = 
\/r^ — 1 cos 9p) and perpendicular [P^/nieC = — 1 sin 9p) to the magnetic field lines 
are, respectively, described as (Harding et al. 2005; Hirotani 2006) 

dn 

dt 
and 



eE\\ - Psc cos 9p, (8) 



dP± „ c dB , , 

— = -i'„sm<), + --P,. (9) 

where 9^ is the pitch angle, P^c represents the radiation drag of the synchrotron and curvature 
radiation, and the second term on the right hand side on equation ([9]) represents the adiabatic 
change along the dipole field line. The radiation drag, Psc? of the synchrotron-curvature 
radiation is described by (Cheng & Zhang 1996), 
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3.3 Model parameters 

The model parameters are the magnetic inclination angle a and the current components 
{ill 32) 1 where ji (or J2) is the current carried by the positrons (or electrons) coming into the 
gap through the inner (or outer ) boundary. We also use the outer boundary as the model 



10 J. Takata and H.-K. Chang 

parameter. As the lower boundary of the gap, we define the last-open field line. For the 
upper boundary of the middle age pulsars, which we discuss in the paper, we use the critical 
field line defined by equation of ([T]). Setting the model parameters, we solve the position 
of the inner boundary, the accelerating electric field, and the current component carried by 
particles created in the gap. 

The created current, and the resultant calculated fiux are determined by the value of 
the injected currents at the inner boundary ji and at the outer boundary j2- In fact, the 
magnitude of the created current is less affected by the current ji injected at the inner 
boundary, but is determined by j2 injected at the outer boundary, because most of pairs 
are created between the ingoing 7-rays and the outward propagating surface X-rays. In 
this paper, we choose the value of the component j2 so that the created pairs produce the 
comparable 7-ray fiux with the observations. 

The created current and the 7-ray fiux do not monotonically increase with the current 
component j2- This is because an increase of the created current strengthens screening of 
the accelerating electric field. With a smaller electric field, the 7-ray photons with a lower 
energy are produced, and therefore a smaller amount of the pairs are created. Namely, the 
increment of the created current is determined under the competition between the increase 
of j2 and the screening of the electric field. For example, if we assume the current j2 = 0.003 
in units of Goldreich- Julian value, we obtain the created current of jg ~ 0.3 for the Geminga 
pulsar with a = 60°. Then if we double the value of j2, we find from the numerical results 
that the increments of the created current and the total luminosity are factors of ~ 1.3 and 
~ 1.2, respectively. We find that the computed 7-ray fiux is not sensitive to the detailed 
values of the current components injected at the boundaries. 

Since the magnetic field must be modified by the rotational and the plasma effects the 
vicinity of the light cylinder (Muslimov & Harding 2005), the last-open field line will be 
also different from the conventional magnetic surface, which is tangent to the light cylinder 
for the vacuum case. For example, Romani (1996) assumed the last-open lines with the 
field lines which have the polar angle 9pc ~ 1.4^/^6'pc,o at the stellar surface, where Opcfl is 
the conventional polar cap angle ([3]) for static dipole field. In the studies of the force-free 
magnetosphere, 9pc ~ 1.36^/^^pc,o was obtained by Contopoulos et al. (1999) and 9pc ~ 
1.27i/2^pc,o by Gruzinov (2005). In the paper, we use 9p^ ~ 1.36i/2^pc,o for the last -open field 
lines. The position of the last-open field lines affects the accelerating electric field, spectrum 
and the pulse profile. For example, if we adopt a vacuum solution, the typical electric field in 
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the gap is described as ~ B^RlSOg^p/cRc with i?* begin the stellar magnetic field (Cheng 
et al. 2000; Hirotani 2006). Then as the last-open field line is located on the magnetic field 
line with a larger polar angle on the stellar surface, we obtain a stronger accelerating electric 
field because (i) the trans-field thickness of the gap, 66 = 6pc — 6c, becomes large and (ii) 
the typical curvature radius of the magnetic field line, on which most 7-rays are emitted, 
becomes small. As a result, a spectrum of the radiation from the gap becomes hard. For 
the pulse profiles, the phase separation between two peaks becomes wide as the last-open 
field shift to inner magnetosphere from the canonical field line, which is tangent to the light 
cylinder (Takata & Chang 2007). In fact, the calculated 7-ray spectrum and pulse profiles 
with the last-open field fines for 6pc ~ 1.36^/2^^^, are more consistent with the observations 
for the Geminga pulsar than those for 6pc ~ ^pc,o (section [5]). 

4 THREE-DIMENSIONAL GEOMETRICAL MODEL 

The observed pulse profile provides an important tool for diagnose the geometry of the 
emission region in the pulsar magnetospheres. In particular, number of the peaks in a single 
rotational period, positions of the peaks in the pulse profile, and morphology of the pulse 
profiles at different energy bands will discriminate the emission models. 

In the framework of the outer gap model (and also the slot gap model), the peaks appear 
in the pulse profiles because the photons emitted various point pileup due to the effect 
of the special relativity correction; the aberration of the emitting direction and the flight 
time (Romani & Yadigaroglu 1995; Dyks & Rudak 2003). To compute the pulse proflle, we 
anticipate that the emission direction coincides with the particle motion. In the observer 
frame, the particle motion along the fleld lines in the north hemisphere is described by 
V = ±VQh + VcoG(j,, where the plus (or minus) sign represents the outgoing (or ingoing) 
particles, vq is the velocity along the magnetic fleld line and is determined by the condition 
that |v| = c, b is the unit vector of the magnetic fleld, for which we adopt a rotating 
dipole fleld in the observer frame, Vco is the co-rotating velocity, and e,^ is the unit vector 
of the azimuthal direction. The emission direction, n = v/c, is interpreted in terms of the 
viewing angle ^ = and the pulse phase $ = — 0e — n ■ Ve/Ric, where 0e and r^ is the 
azimuthal emission direction and the radial distance to the emission point, respectively. We 
use a rotating magnetic dipole fleld in the observer frame to compute the pulse proflle. 

The distribution of emissivity affects the pulse proflles. For the distribution of the emis- 
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sivity in the meridional plane, we use the results of the two-dimensional electrodynamic 
model, which solves the trans-field distribution of the emissivity in the magnetic merid- 
ional plane (see section [5.1.3p . For the structure of the azimuthal direction, it is difficult 
to model without solving a three-dimensional dynamics. For the first order approximation, 
however, we may apply the function form of the accelerating electric field for the vacuum 
solution. In the vacuum solution, the typical value of the electric field along a field line 
is determined by oc Rn{4>)~^ with Rn being the curvature radius at the null charge 

surface and </> the magnetic azimuthal angle of the magnetic field on the stellar surface 
measured from the magnetic meridional plane (Cheng el al. 1986a; Cheng et al. 2000). In 
such a case, the number of the curvature photons emitted per second by a particle and the 
typical energy of the emitted photons are proportional to Nc{(f)) oc Rn^^^ and Ey oc Rn^^^-, 
respectively, if we assume the force balance between the particle acceleration and the radia- 
tion drag. For example, we obtain Rn ~ 0.24i?ic for the magnetic meridional plane (</> = 0) 
and Rn ~ 0.72Ric for = 90°, if we use 9pc ~ 1.36^/^6'pc,o for the last-open field lines with 
the inclination angle a = 60°. In this paper, therefore, when we compute the pulse profile, 
we assume that ratio of the photon number emitted on the field lines of = 0° and of 
(f) = 90° is [iV(90°)/iV(0°)] = (0.72/0.24)-3/4 ~ 0.42. We do not consider the emissions from 
the magnetic field lines on which the radial distance to the null charge point is larger than 
r = Ric- 

5 RESULTS 

5.1 The Geminga pulsar 

In this section, we apply the model to the Geminga pulsar and constrain the magnetic 
inclination and viewing angles. First, we will show structure of the outer gap that occupies 
most of the region between the last-open line and the critical field line (section 15.1. II) . Then 
we will show that the calculated spectra are consistent with the observations from X-ray 
through 7-ray bands (section [5.1.2p . We use the properties of the observed 7-ray pulse profile 
to constrain the magnetic inclination and viewing angles of the observer (section [5.1.3l) . For 
the Geminga pulsar, the observed pulse profile in 7-ray bands shows double peak structure, 
and the phase separation of the two peaks is about 0.5 phase. Finally, we further constrain 
the magnetic inclination angle by discussing visibility of the radio beam in section 15.21 
in which we will assume that the radio beam from the Geminga pulsar is oriented in a 
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Figure 2. Left; The geometry of the outer gap accelerator. The soHd Une is the boundary of the gap. The dotted lines show 
the critical line, and the dashed-dotted line is the null surface of the Goldreich-Julian charge density. Right; The trans-field 
structure of the current (solid) and the electric potential drop (dashed line) along the magnetic field lines as a function of 
the polar angle of the field lines on the stellar surface. We note that our static approximation will not be good treatment to 
describe the electric structure beyond the light cylinder. But we expect that the gap thickness and the calculated spectrum 
will not be affected much the electric structure beyond light cylinder, because most 7-ray photons are radiated from the lower 
part outer gap inside the light cylinder. The results are for 9pc = l-SG'^/^SpcO- 




Figure 3. The distribution of the accelerating electric field along the magnetic field line, on which the strength of electric field 
attains the maximum value. The solid and the dashed lines are results from the last-open field lines for 9pc = 1.36^^^9pcfi and 
9pc = l■36^''■^^?pc,0I respectively. 



different direction from the line of siglit. For tlie Geminga pulsar, we find that the calculated 
spectrum with the stellar dipole momentum inferred from the rotating dipole model is too 
soft compared with the observation. For the present calculations, we used double of the 
dipole momentum inferred from the rotating dipole model. 



14 J. Takata and H.-K. Chang 
5.1.1 Outer gap structure 

Figure [2] summarizes the outer gap structure of the Geminga pulsar for the inchnation angle 
that a = 60°, and for ji 2 = 3-10"^, that is 0.3 % of the Goldreich- Julian currents are injected 
at both inner and outer boundaries. For the last-open field line, we use the magnetic field 
line that has the polar angle 9pc ~ 1.36^^'^6pcfl, as we discussed in section 13^31 

In the left panel in Figure [2], the thick solid line shows the boundary of the obtained 
gap structure. The dotted and dashed-dotted line represent the critical field line and the 
null charge surface of the Goldreich- Julian charge density, respectively. As the computed 
boundary of the gap shows, we can find the solution that the outer gap can occupy whole 
region between the last-open field line and the critical field line, as we discussed in section |2j 

With the right panel of Figure [2], we display the current density carried by the created 
pairs (solid line) and the total electric potential drop in the gap (dashed line), that is 
potential drop between the inner and outer boundaries, as a function of the polar angle, 6, 
of the magnetic field lines on the stellar surface. In the abscissa, the polar angle satisfying 
the condition that {6pc — 0)/ 59 gap = is corresponding to the last-open field line, and the 
angle of {6pc — 9)/56gap = 1 is corresponding to the upper boundary of the gap. From the 
solid line, we can see that the maximum current density of about 30 % of the Goldreich- 
Julian value runs through along the magnetic field line that penetrates the gap at the height 
of about 40% of the gap thickness measured from the last-open field line. From the dashed 
line, we also find that the largest electronic potential drop oi V ^ 2 x 10^^ Volt appears 
on the magnetic field line penetrating the gap at the height of about 20% of the thickness. 
This indicates that the strength of the accelerating electric field becomes maximum on the 
magnetic field line penetrating the gap at the height of about 20% of the thickness. Because 
the emitting power of each magnetic flux tube is proportional to the product of the current 
density and the potential drop at each field line, we can read that most of emissions come 
from the height between 20 % and 40 % of the gap thickness (see also Figure [5]). 

Figure [3] shows the distribution of the electric field along the magnetic field line, on which 
the strength of the electric field attains the maximum value in the outer gap. We compare 
the results for 6pc = l-36^/^6'pc,o (solid line) and 9pc = 9pcfi (dashed line). As the Figured 
shows, a larger accelerating electric field is obtained if we assume a larger polar cap angle. 
The typical strength of the accelerating electric field is proportional to E\\ oc [SOgap)"^ /Rc 
(Cheng et al 2000; Hirotani 2006). In the present case, the curvature radius on the field 
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Energy [GeV] 

Figure 4. The curvature-synchrotron spectrum for the Geminga pulsar in X-ray and 7-ray bands. The solid line represents 
the total emissions which includes the outward (dashed line) and the inward (dotted lines) emissions. The filled circles and the 
dashed-dotted line represent the observed phasc-avcragcd spectrum (after Kargaltsev et al. 2005). 

line on which the electric field attains the maximum value is typically Rc ~ O.SRic for 
9pc = 1.36^/^^pc,o and ~ l-2Ric for 9pc = Opcfl- This means that the ratio of strengths 
of the accelerating electric field for two cases is E\\\q^^^i^^qi/2q^^ ^/ E\\\Q^^=Q^^^^ ~ 2, which 
describes the difference of the strengths of the electric field seen in Figure [31 The typical 
energy of the photons emitted by the curvature radiation is proportional to Ec oc E'^^^R].^'^. 
The ratio of the typical energies becomes Ec\q^^=i^^qi/2q^^ Ec\epc=ep^,o ~ 1-4. As a result, the 
calculated spectrum for 9pc = 1.36^/^^pc,o is harder than that for 9pc = 9pcfi. 

5.1.2 Spectrum 

Figure H] compares the calculated spectrum of the curvature and synchrotron processes with 
the observations in X-ray through 7-ray bands. The result is for the last-open field lines 
with the polar angle 9pc = 1.36^/^6'pc,0; and the model parameters are the same with that of 
Figure El In Figure HI the filled circles and the dashed-dotted line show the observed phase- 
averaged spectrum. The dashed line and the dotted line represent the emissions of outgoing 
positrons and of the ingoing electrons, respectively. The solid line represents the total emis- 
sions. For each line, the components of the curvature and the synchrotron radiations are 
combined. In fact, the curvature radiation dominates in the spectrum above 100 MeV, while 
the synchrotron radiation dominates below 10 MeV. From Figure HI we see that the total 
spectrum (solid line) is consistent with the observations in X-ray through 7-ray bands. 
If we use the canonical last-open field line with the polar angle 9pc = 9pc^Q, the spectral 
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Figure 5. The trans-field distribution of the emissivity of the 7-ray emissions. Each hne represents the spectrum from the 
height below 20 % of the thickness measured from the last-open field line, (solid line), between 20 % and 40 % (dashed line), 
between 40 % and 60 % (dotted line), and above 60 % (dashed-dottcd line). 



cut off is lactated at the energy factor of 1/1.4 smaller tlian tliat seen in tlie solid line of 
Figure HJ as we discussed in the last paragraph of section I5.1.1[ Therefore, the calculated 



the canonical last-open field line. 

From Figure HJ we can see that the curvature radiation of outgoing positrons (dashed 
line) dominate in the total emissions above 100 MeV. This is because most of the pairs 
are created the vicinity of the inner boundary by the inward propagating 7-rays and the 
outward propagating surface X-rays. The created electrons and positrons are accelerated 
inwardly and outwardly by the electric field, respectively. The positrons created around 
the inner boundary are accelerated by almost whole potential drop in the gap, while the 
created electrons are accelerated by the potential drop only between the inner boundary 
and the point where they were created. Therefore, the total radiation power is larger for the 
outward emissions than the inward emissions. Below 10 MeV, on the other hand, we find 
that the synchrotron radiation of the ingoing electrons dominate in the calculated spectrum. 
This is because the synchrotron emissions with strong magnetic field near the stellar surface 
contribute to the spectrum in X-ray bands. 



spectrum for the last-open field line with 6, 



pc 



1.36^/^6'pc,o is more consistent with that for 



5.1.3 '-) -ray pulse profiles 

We discuss the 7-ray pulse profiles, and constrain the inclination angle and the viewing angle 
using the three-dimensional geometrical model described in section |4l As Figure H] shows, the 
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Figure 6. Expected pulse profiles in the 7-ray bands for the Geminga pulsar. The upper left, the upper right, the lower left, 
and the lower right panels show the pulse profiles for the inclination angle of q = 20°, 40°, 60° and 80°. The solid and dashed 
lines in each pane represent the pulse profiles of the viewing angle 5 = 85° (or 95°) and 70° (or 110°), respectively. The results 
are for Ope = l.SO^/^g^^ ^ 



outward emissions dominate in the total spectrum above 100 MeV. In this section, therefore, 
we consider only the outward emissions to compute a pulse profile. 

For the trans-field distribution of the emissivity in the poloidal plane, we use the results 
of the two-dimensional electrodynamic study. In Figure O we can read the trans- field distri- 
bution of the emissions above 10 MeV bands. For example, the dashed line represents the 
spectra of the emissions from the height between 20 % and 40 % of the thickness measured 
from the last-open field line. By comparing the magnitude of the fluxes between the solid 
and dashed lines in Figure [5l we may put an emissivity at the height below 20% of the gap 
about one order smaller than that at the height between 20% and 40% of the gap. 

Figure [6] summarizes the pulse profiles for the magnetic inclination angle, a = 20° (upper 
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Figure 7. The pulse profile of the Geminga pulsar with the inclination angle a = 60° and the viewing angle ^ = 85° (or 105°). 
The solid and dashed lines represent the pulse profiles for the outward and inward emissions, respectively. Above 100 MeV 
bands, only outward emissions contribute to the spectrum (Figure |4]l so that the pulse profile (solid line) becomes double peak 
structure. For around 10 MeV bands, both the outward and inward emissions contribute to the pulsar profile, and the pulse 
profile has the four peaks in a single period. For X-ray bands, the inward emissions dominate in the pulse profile, which has 
two peaks (dashed line) in a single period. 

left panel), a = 40° (upper right panel), a = 60° (lower left panel) and a = 80° (lower right 
panel). We calculate the emissions from the whole region between the last-open field lines 
and the critical field lines. The results are for the last-open field lines with the polar angle 
Opc = If we used the canonical last-open field lines with 9pc = 9pc to computed 

pulse profile, we found that no geometries (inclination angle and viewing angle) can produce 
the observed phase separation between the two peaks with the 7-ray spectrum; we obtained 
narrower phase separation {6(f) < 0.4) than the observation 6(f) ~ 0.5. 

In Figure [6], the solid and dashed lines are results for the viewing angle of = 85° (or 
95°) and ^ = 70° (or ^ = 110°). As the same with the conventional outer gap model, the first 
peak in the pulse profile is made by the emissions from leading part of the gap according to 
the rotational direction, while the second peak is made by the emissions from trailing part of 
the gap. We can see from Figure [6] that the morphology of the pulse profiles depends on the 
inclination and viewing angles. By comparing the observed pulse profile, therefore, we can 
constrain the inclination and viewing angles. First we can rule out nearly aligned rotator, 
as represented by the pulse profile for a = 20°. With such a small inclination angle, a sharp 
second peak does not appear in the pulse profile for any viewing angles of the observer. 
Secondary, we would rule out the viewing angle far away from C, = 90° as demonstrated by 
the pulse profile (dashed line) with C, = 70° in Figure [61 In such a viewing angle, the pulse 
profile has a single broad peak or the two peaks with a narrower separation of the phases 
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compared with the observations (~ 0.5 phase). Therefore, a larger inchnation angle and the 
viewing angle that ^ ~ 90°, say 80° < ^ < 100°, are preferred to explain the observed pulse 
profile of 7-ray bands. 

5.1.4 Pulse profiles below 10 MeV 

In this section, we discuss the expected pulse profile in soft 7-ray and the X-ray bands. 
As the calculated spectrum in Figure H] shows, the inward emissions from the inward parti- 
cles contribute to the total spectrum below 10 MeV. Specifically, both the inward and the 
outward emissions contribute to the total emissions around 10 MeV, while the only inward 
emissions dominate into the total emission in X-ray bands. Therefore, we expect that the 
pulse profile below 10 MeV is not simply an extension of that of the 7-ray bands above 
100 MeV. Figure [7| shows the pulse profile for the outward (solid) and the inward (dashed) 
emissions for a = 60° and ^ = 85° (or 95°). From the dashed line in Figure [3, we find that 
the pulse profile of the inward emissions has also double peak structure. For the inward 
emissions, the first peak is made by the emissions from the trailing part of the gap and the 
second peak is made by the emissions of the leading part of the gap. 

For the calculated spectrum around 10 MeV in Figure HI both the outward and inward 
emissions contribute to the spectrum. The pulse profiles in 10 MeV bands will be made by 
both the outward (solid line in Figure [7j) and inward (dashed line in Figure [7j) emissions. 
Therefore, we expect that the pulse profile in soft 7-ray bands will have four peaks in a 
single period. 

In X-ray bands, the inward emissions dominate in the total emissions (Figure H]). There- 
fore we expect the double peak structure of the pulse profile in the X-ray bands for the 
Geminga pulsar. By comparing between the predicted pulse profiles in 7-ray (solid line) 
above 10 MeV bands and X-ray bands (dashed line), the pulse peaks in X-ray band will 
appear in the leading phases of the pulses of the 7-ray bands. 

With the present results, the observed differences of the pulse morphologies in the dif- 
ferent wavelength are strong tool to diagnose the emission processes in the pulsar mag- 
netosphere. For the Vela pulsar, the outer gap model can explain the morphology differ- 
ences between the pulse profiles in RXTE and EGRET bands (Takata et al. 2008). For 
the Geminga pulsar, the spectrum and the pulse profiles observed by ROAST and ASCA 
indicated that the non-thermal X-rays and 7-rays are produced by a distinct mechanism 
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(Jackson et al. 2002). Kargaltsev et al. (2005) analyzed X-ray spectrum and pulsations ob- 
served by XMM — Newton. They found that the pulse profile of the non-thermal X-rays in 
2-8 keV bands has two peaks in a single period, and that the spectrum and pulse profile are 
not simply an extension of that of the EGRET bands. The observed pulse morphologies for 
the Geminga pulsars are different in the different wave length. As we discussed above, the 
present model also predicts the double peak structure in X-ray pulse profile and the energy 
dependent pulse profiles from X-ray through 7-ray bands. The results of the present model 
are generally consistent with the observed pulse profiles in X-ray and 7-ray bands. For the 
phase separation between two peaks in X-ray bands, however, the observations in 2-8 keV 
bands imply about 0.35 (0.65) phase (Kargaltsev et al. 2005), while the present model pre- 
dicts about 0.5 phase (dashed line in Figure [7]). So, it may be required an examination of 
the X-ray pulse profile with a more realistic magnetic field configuration. 

5.2 Radio "quiet" observations 

The Geminga pulsar is the only known radio quiet (or weak emissions) 7-ray pulsar. Based 
on the following reasons, we assume that the radio beam of the Geminga pulsar is oriented 
in a different direction from the line of sight. First, there are many radio pulsars around 
the Geminga pulsar on P — P diagram. Second, Lin & Chang (2008) searched candidate 
periods in unidentified EGRET sources for the next Geminga. They found that most of the 
candidate timing properties belong to that of Vela-like pulsars. If their candidate periods are 
true, their results would suggest that the major cause for a radio quiet is likely geometrical 
rather than intrinsic. 

The radio emission model with the assumption that the radio "quiet" for the Geminga 
pulsar is caused by geometry also provides a tool to study the magnetic inclination angle and 
the viewing angle with the 7-ray emission model. For the emission height of the outer cone 
of the radio emissions, we apply the empirical model developed by Kijak & Gil (2003). They 
found a relation among the emission height, the rotation period, the period time derivative 
and the frequency of the radio wave. Assuming the rim of the radio beam are emitted vicinity 
of the critical magnetic field line (defined by 6'c), their relation is described as 

r ~ mR.{epc/e.)'p'-\p/iQ-'h,-'r'v^^\ (14) 

where i?* is the stellar radius and v is the frequency of the radio wave. 

Figure [8] shows photon mapping of the rim of the radio beam emitted at the height 
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Figure 8. Photon mapping of the outer rim of the radio beam emitted at the height described by equation of II14I I for different 
inclination angles a = 20° (upper left), 40° (upper right), 60° (lower left) and 80° (lower right). The dotted line represent 
viewing angle of a = 85°. 



described by equation of (fT4|) for different inclination angles. We used the rotating dipole 
field as the magnetic field of the magnetosphere. In Figure [8], the solid and dashed lines 
represent the emissions for the frequencies of z/ = 100 MHz and 400 MHz, respectively. For 
the references the dotted line represents line of sight for the viewing angle of ^ = 85°. 

For a small inclination angle represented by a = 20° (upper left panel in Figure [8]), 
an observer with S, ~ 90°, which is preferred to explain observed 7-ray pulse profile, does 
not measure the radio emission. This is because the magnetic pole points in a completely 
different direction from the line of sight for ^ ~ 90°. For largely inclined pulsar such as 
a = 80° (lower right panel), on the other hand, it is highly possible to be measured the 
radio emissions by the observer with ^ ~ 90°, because the magnetic pole points close to the 
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Table 1. Timing properties of the middle-aged pulsars. P is the rotational period, P is the time derivative of the rotational 
period, B is the stellar magnetic field inferred from the dipole field model, r is the pulse age, d is the distance, E/d? is the 
spin-down luminosity measured on Earth, a is the magnetic inclination angle, and ^ is the viewing angle. For the magnetic 
inclination angle a and the viewing angle ^, we refer the present paper for the Geminga pulsar, Rankin (1993) for PSRs 
B0355+54, B0740-28 and B1449-64, and Everett and Weisberg (2001) for PSR B1929+10. In the last column, we show the 
predicted detectability of pulsed 7-ray emissions by Fermi telescope with 5o- in 50 hours. 

line of sight. For the Geminga pulsar, therefore, a smaller inclination angle is preferred to 
explain the radio "quiet" observations with the 7-ray detections. 

Using the observed pulse profiles in the 7-ray observations discussed in section 15.1.31 
we safely eliminated the possibility of the nearly aligned rotator such as a = 20° because 
predicted single peak structure (upper left panel in Figure [6]) is not consistent with the ob- 
served double peak structure of the 7-ray pulse. Combing those results obtained by analysis 
in the 7-ray and the radio "quiet" observations, we conclude that the Geminga pulsar has a 
moderate magnetic inclination angle of a ~ 50°, and we observe the Geminga pulsar with a 
viewing angle of ~ 90°. 



5.3 Another middle-aged pulsars 

We apply the present outer gap model to another middle-aged pulsars. Using ATNF Pulsar 



Catalogue (Manchester et al. 2005, see http://www.atnf.csiro.au/research/pulsar/psrcat), 
we chose four middle-aged pulsars, PSRs B0355+54, B0740-28, B1449-64, and B1929+10, 
to discuss the detectability of the 7-ray emissions from the outer gap by Fermi telescope. 
These four pulsars were chosen with the following conditions. If the pulsars are at a large 
distance from the Earth, Fermi telescope could not measure the 7-ray emissions from the 
pulsars. A rough idea of the 7-ray flux on the Earth is described by r]Lsp/(P, where 77 is 
efficiency, L^p is the spin-down luminosity, and d is the distance to the pulsar. Our four 
pulsars have the spin-down luminosity so that rjLsp/d'^ with 77 = 0.01 of the Geminga pulsar 
is larger than 10~^^ erg/cm^s, which is close to the lower limit of Fermi telescope detection 
with 5(7 in 50 hours. 

Second, there is the possibility that we can not measure the 7-ray emissions, because the 
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Figure 9. The 7-ray spectrum (left) and pulse profile (right) for PSR B0355+54. For a = 51°, the observe with ^ ~ 55° 
measures only the outward emissions. 

7-ray beam is oriented in a different direction from tlie line of sight. Because tlie geometry 
of the 7-ray beams depends on the emission models, the detectability of the 7-ray emissions 
depends on the models. If we know the viewing angle and the magnetic inclination angle 
of the pulsar, we can test the models with Fermi observations. On these ground, we chose 
those four pulsars because the magnetic inclination and viewing angles are inferred from the 
radio polarization data (Rankin 1993; Everett & Weisberg 2001). In Table 1, we summarize 
the timing properties of the four pulsars and the Geminga pulsar. The represented results 
in this section were computed with (i) the double of the dipole inferred from the rotating 
dipole model and with (ii) the last-open field lines with the polar angle 6pc = 1.36^/^6'pc,o- 



5.3.1 PSR BO355+54 

PSR B0355+54 has a rotational period of P = 0.1564 s (Manchester et al. 1972), which is 
the fastest in the four pulsars we chose. The distance from the Earth is ~ 1.1 kpc (Hobbs 
et al. 2004). PSR B0355-h54 is known to emit the X-ray (McGowan et al. 2006) and probably 
TeV 7-rays (Senecha et al. 1995) from its pulsar wind nebula. The polarization data on the 
radio wave can be fitted by the rotational vector model; the magnetic inclination and the 
viewing angles are inferred as a ~ 51° and ^ ~ 55°, respectively (Rankin 1993). 

Figure [9] shows the expected spectrum (left panel) and the pulse profile (right panel) in 
7-ray bands for PSR B0355+54. For the current components injected at the inner and outer 
boundaries, we used ji = j2 ~ 0.01. For the inclination angle a = 51°, the observer with the 
viewing angle ^ = 55° can not measure the inward emissions from the outer gap. The model 
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predicts the emission flux around 1 GeV are about F ~ 10~^^ erg/s, which will be able to 
be measured by Fermi telescope. With the viewing angle ^ ~ 55°, the expected pulse profile 
has a broad peak in a single rotational period, as right panel in Figure [9] shows. 

5.3.2 PSR B0740-28 

PSR B0740-28 is younger and faster rotating pulsar than the Geminga pulsar (Bonsignori- 
Facondi et al 1973). The distance to the pulsar is d = 1.89 kpc (Hobbs et al. 2004) and the 
magnitude of L^p/d"^ is the largest in the four pulsars. The pulsar wind powered by PSR 
B0740-28 was discovered in optical bands (Jones et al. 2002). The magnetic inclination angle 
a ~ 37° and the viewing angle ^ ~ 27° are inferred with the polarization characteristics of 
the radio wave (Rankin 1993). 

We will not expect the detection of the 7-ray beam by Fermi telescope from PSR B0740- 
28, because the 7-ray beam will be oriented in a different direction from the line of sight. 
Because the viewing angle with ^ ~ 27° points between the rotational and magnetic axes 
in the magnetic meridional plane, the observer measures mainly emissions inside the null 
charge surface and on the magnetic field lines toward rotational axis. However, the outer gap 
does not predict a strong emissions from that region, but predicts a strong emission beyond 
the null charge surface and on the magnetic field lines curved away from the rotational axis. 
The present outer gap model expects that the strong emissions from the outer gap of PSR 
B0740-28 is oriented in a different direction from the line of sight. 

For the slot gap model (Harding et al. 2008), on the other hand, a detection of 7-ray 
emissions from PSR B0740-28 may be expected, because the model assumes the emissions 
from the polar cap region to the vicinity of the light cylinder above the whole last-open 
field lines (Dyks & Rudak 2003). Therefore, PSR B0740-28 will be a good candidate to 
discriminate the outer gap model and the slot-gap model with Fermi telescope. 

5.3.3 PSR B 1449-64 

PSR B1449-64 has the rotational period P = 0.1785 and the period time derivative P = 
2.746 ■ 10~^^ s/s (Large et al. 1969). The distance to the pulsar is d = 1.84 kpc (Siegman et 
al. 1993). The magnetic inclination of a ~ 56° and the viewing angles of ^ ~ 63° are used 
to fit the polarization data of the radio wave with the rotating vector model (Rankin 1993). 
Figure [10] shows the expected spectrum (left panel) and the pulse profile (right panel) in 
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Figure 10. The 7-ray spectrum (left) and pulse profile (right) for PSR B1449-64.In the left panel, the solid, dashed and dotted 
lines represent the spectra of the total, outward and inward emissions. In the right panel, the solid and dashed show the pulse 
profile of the outward and inward emissions, respectively. 
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Figure 11. The 7-ray spectrum (left) and pulse profile (right) for PSR B1929+10.The lines correspond to same case as 
Figure [To] 



7-ray bands for PSR B1449-64. The expected pulse profile shows the double peak structure. 
We adopt ji = j2 = 0.02 for the current components injected at the inner and the outer 
boundaries. 



5.3.4 PSR B 1929+ 10 

PSR B1929+10 is known as one of the closet pulsar with the distance d ~ 0.36 kpc (Hobbs 
et al. 2004). PSR B1929+10 is the oldest pulsar r ~ 3 ■ 10^ yrs in the pulsars we chose with 
its rotation period P = 0.2265 s and the period time derivative P = 1.157 ■ 10~^^ s/s. The 
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X-ray emissions from both the pulsar and the pulsar wind nebula have been reported (Becker 
et al. 2006; Hui and Becker 2008). The inferred magnetic inclination and the viewing angle 
of the observer from the radio polarization are a ~ 35° and ~ 60°, respectively (Everett 
and Weisberg 2001). We use ji = j2 = 0.1 for the currents injected the inner and outer 
boundaries. 

Because of the closeness to the Earth of PSR B1920+10, the expect 7-ray flux can signif- 
icantly exceed the lower limit of Fermi telescope detection with 5a in 50 hours (Figure [TTi) . 
although the spin-down power and therefore the 7-ray luminosity are the smallest in the 
four pulsars. For the pulse profile, we expect single peak in a single rotation period. 

6 SUMMARY 

We summarize the results of our work as follows. 

1. We have studied the non-thermal emissions of the Geminga pulsar with the outer gap 
accelerator model. We solved the accelerating electric field with the pair creation and the 
radiation processes in the magnetic meridional plane. We demonstrated that the calculated 
spectrum of the curvature and the synchrotron emissions is consistent with the observation 
in X-ray through 7-ray bands. The outward curvature radiation dominates in the spectrum 
of 7-ray bands, while the inward synchrotron emissions dominate in X-ray bands. 

2. We demonstrated that the 7-ray pulse above 100 MeV has two peaks, which are made 
by the outward emissions. Around 10 MeV bands, we expect four peaks in a single period 
with both the outward and inward emissions. In X-ray bands, the pulse profile has two 
peaks made by the inward emissions. The pulse phase in X-ray bands are not in phase with 
that in the 7-ray bands. The model predicts the energy dependent pulse morphology in X- 
ray through 7-ray bands, which is also indicated by the observations. (Jackson et al. 2002; 
Kargaltsev et al. 2005). 

3. We discussed the inclination angle and the viewing angle for the Geminga pulsar with 
the pulse morphology. The observed double peak structure in EGRET bands rules out nearly 
aligned rotator. Assuming that the radio beam from the Geminga pulsar is oriented in a 
different direction from the line of sight, a greatly inclined magnetic axis is also ruled out. 
The observed phase separation of two peaks in 7-ray pulse profiles is explained by a viewing 
angle C, ~ 90°. We conclude that the Geminga pulsar has a inclination angle that a ~ 50° 
and a viewing angle that C, ~ 90°. 



Non-thermal emissions from middle-aged pulsars 27 
4. We applied our method to another four middle-aged radio pulsars, whose spin-down 
power and distant from the Earth expect the potential of the detection of the 7-ray emissions 
from those pulsars by Fermi telescope. Applying the inclination angle and the viewing angle 
inferred from the radio polarization characteristics, we discussed the visibility of the 7-ray 
emissions from the outer gap. We predict that 7-ray emissions from PSRs B0355-I-54, B1449- 
64 and B1929+10 are detectable by Fermi telescope. For PSR B0740-28, the 7-ray beam 
from the outer gap will be oriented in a different direction from the viewing angle. The 
observations of the PSR B0740-28 by Fermi telescope will probably be useful for studying 
the emission site of not only the 7-rays but also the radio wave in the pulsar magnetospheres. 

ACKNOWLEDGMENTS 

The authors appreciate fruitful discussion with K.S. Cheng, K.Hirotani, K.S. Cheng, S. 
Shibata, and R.Taam. Authors also thank anonymous referee for his/her insightful comments 
on the manuscript. This work was supported by the Theoretical Institute for Advanced 
Research in Astrophysics (TIARA) operated under Academia Sinica and National Science 
Council Excellence Projects program in Taiwan administered through grant number NSC 
96-2752-M-007-002-PAE. 

REFERENCES 

Arons J., 1983, ApJ, 266, 215 

Becker W., et al., 2006, ApJ, 645, 1421 

Bonsignori-Facondi S.R., Salter C.J. & Sutton J.M., 1973, A&A, 27, 67 

Chang H.-K., Boggs S., Chang Y.-H. for the NCT collaboration, 2007, AdSpR, 40, 1281 

Cheng K.S., Ho C. & Ruderman M. 1986a, ApJ, 300, 500 

Cheng K.S., Ho C. & Ruderman M. 1986b, ApJ, 300, 522 

Cheng K.S., & Zhang L. 1996, ApJ, 463, 271 

Cheng K.S., Ruderman M. & Zhang L. 2000, ApJ, 537, 964 

Contopoulos I., Kazanas D. & Fendt C. 1999, ApJ, 511, 351 

Daugherty J.K. & Harding, A.K., 1996, ApJ, 458, 278 

Dean et al., 2008, Sci, 321, 1183 

Dyks J. & Rudak B., 2003, ApJ, 598, 1201 

Dyks J., Rudak B. & Harding A.K., 2004, ApJ, 607, 939 



28 J. Takata and H.-K. Chang 
Everett J.E.& Weisberg J.M., 2001, ApJ, 553, 341 

Fierro J.M., Michelson P.F., Nolan P.L. & Thompson D.J., 1998, ApJ, 494, 734 
Grizov A., 2005, Phys. Rev. Lett. 94, 021101 

Harding A.K., Stern J.V., Dyks J. & Frackowiak M., 2008, ApJ, 680, 1378 
Harding A.K., Strickman M.S., Gwinn C, Dodson R., Moffet D. & McCulloch P., 2002, 
ApJ, 576, 376 

Harding A.K., Usov V.V., Muslimov A.G., 2005, ApJ, 622, 531 
Hirotani K., 2006, ApJ, 652, 1475 
Hirotani K., 2007, ApJ, 662, 1173 

Hobbs G., Lyne A.G., Kramer M., Martin C.E. & Jordan C, 2004, MNRAS, 353, 1311 
Hui C.Y. & Becker W., 2008, A&A, 486, 485 

Jackson M.S., Halpern J. P., Gotthelf E.V. & Mattox J.R., 2002, ApJ, 578, 935 
Jones D.H., Stappers B.W.& Gaensler B.M., 2002, A&A, 389L, 1 

Kanbach G., Slowikoska A., Kellner S. & Steinle H., 2005,2005, in Bulik T., Rudak B., 
Madejski G., eds, AIP Conf. Proc. Vol. 801, Astrophysical Sources of High Energy Particles 
and Radiation. American Institute of Physics, New York, p. 306 Proceeding, 801, 306 

Kargaltsev O.Y., Pavlov G.G., Zavlin V.E. & Romani R.W., 2005, ApJ, 625, 307 

Kataoka J., et al., 2005, Proc. SPIE, 5898, 133 

Kijak J.& Gil J. 2003, A&A, 397, 969 

Large M.I., Vaughan A.E. & Wielebinski R., 1969, Natur, 223, 1249 
Lin L.C.-C. & Chang H.-K., 2008, MNRAS, 387, 729 
Manchester R.N., Hobbs G.B., Teoh A., Hobbs M., 2005, AJ, 129, 1993 
Manchester R.N., Taylor J.H.& Huguenin G.R., 1972, Nat., 240, 74 

McGowan K.E., Vestrand W.T., Kennea J.A., Zane S., Cropper M.& Cordova F.A., 2006, 
ApJ, 647, 1300 

Mestel L., 1999, Stellar Magnetism, International Series of Monographs of 

Muslimov A.G. & Harding A.K. 2003, ApJ, 588, 430 

Muslimov A.G. & Harding A.K. 2004, ApJ, 606, 1143 

Muslimov A.G. & Harding A.K. 2005, ApJ, 630, 454 

Rankin J.M., 1993, ApJS, 84, 145 

Romani R.W. & Yadigaroglu I.-A. 1995, ApJ, 438, 314 

Romani R.W. 1996, ApJ, 470, 469 

Ruderman M.A. & Sutherland P.G., 1975, ApJ, 196, 51 



Non-thermal emissions from middle-aged pulsars 29 

Senecha V.K., Bhat C.L., Rawat H.S., Rannot R.C, Sapru M.L., Kaul R.K. & Tickoo 
A.K., 1995, A&A, 302, 133 
Shibata S., 1995, MNRAS, 276, 537 

Siegman B.C., Manchester R.N.& Durdin J.M., 1993, MNRAS, 262, 449 

Tang RS. Anisia, Takata J., Jia, J.J., Cheng K.S., 2008, ApJ, 676, 562 

Takata J., Shibata S. & Hirotani K. 2004, MNRAS, 354, 1120 

Takata J., Shibata S., Hirotani K. & Chang H.-K. 2006, MNRAS, 366, 1310 

Takata J. & Chang H.-K., 2007 ApJ, 670, 677 

Takata J., Chang H.-K. & Cheng K.S., 2007 ApJ, 656, 1044 

Takata J., Chang H.-K. & Shibata S., 2008, MNRAS, 386, 748 

Thompson D.J., 2004, in Cheng K.S., Romero G.E., eds. Cosmic Gamma Ray Sources. 

Dordrecht, Kluwer, p. 149 
Vats H.O., Singal A. K., Deshpande M.R., Iyer K.N., Oza R., Shah C.R. & Doshi S., 1999, 

MNRAS, 302, L65 



